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Optical activity is ubiquitous across natural and artificial media and is conventionally understood 
in terms of scattering from electric and magnetic moments. Here we demonstrate experimentally 
and confirm numerically a type of optical activity that cannot be attributed to electric and magnetic 
multipoles. We show that our observations can only be accounted for by the inclusion of the toroidal 
dipole moment, the first term of the recently established peculiar family of toroidal multipoles. 


Chirality, i.e. the property by which a pattern can¬ 
not be superimposed with its mirror image [1], is a ba¬ 
sic characteristic of structured matter at all scales, from 
the cosmological to the molecular. It enters fundamen¬ 
tal questions, such as the homochirality observed in life 
[2] and the parity violation within the Standard Model 
of physics ESI. An ubiquitous manifestation of chiral¬ 
ity is optical activity, an optical effect characterized by 
the rotation of the polarization plane of linearly polarized 
light propagating through a chiral medium, and preferen¬ 
tial absorption of circularly polarized light of a particular 
handedness (circular dichroism) if the medium is dissipa¬ 
tive. Observed more than 150 years ago by Louis Pasteur 
on solutions of tartaric acid [5] , optical activity occupies 
a central position in the diagnostic methodology of many 
scientific disciplines [6], where it provides means of ob¬ 
taining information about the microscopic structure and 
electromagnetic excitations of a medium from its far-held 
scattering properties. An examination of optical activ¬ 
ity is invaluable, for example, in the biosciences where 
it enables deduction of the conformation of proteins in 
a dynamic fashion [7], or in the study of magnetic phe¬ 
nomena in strongly correlated electronic systems [8]. 

Optical activity is typically described within the dipole 
approximation of the dynamic multipole expansion [6], 
where the interaction of an incident wave with collinear 
electric and magnetic dipoles induced in the medium 
affects the polarization state of the wave upon trans¬ 
mission. Certain combinations of electric dipole and 
other higher-order multipoles (most notably electrical 
quadrupole) can also contribute to optical activity jOlfTO]. 
although such contributions are usually considered neg¬ 
ligible, especially in isotropic media [llllig. It was sug¬ 
gested recently m that optical activity could arise even 
in the absence of electric dipolar response due to the 
excitation of toroidal multipoles, an elusive part of the 


dynamic multipole response [Ml [15]. In particular, it was 
anticipated that toroidal dipole should contribute to op¬ 
tical activity in the same manner as electric dipole m 
since the two have the same angular momentum and par¬ 
ity properties mils]. 

In this Letter we observe experimentally and confirm 
numerically that the toroidal dipolar excitation can in 
fact provide the dominant contribution to optical ac¬ 
tivity, assuming the role of the electric dipole. Using 
the metamaterial approach, we have been able to en¬ 
hance in a resonant fashion the microscopic toroidal and 
electric quadrupole responses to circularly polarized ra¬ 
diation and show that this previously unexplored and 
rather exotic combination of dynamically induced mul¬ 
tipoles leads to strong circular dichroism. Our results 
challenge the common interpretation of optical activity 
based on the electric dipole-magnetic dipole approxima¬ 
tion and enhance its capabilities as a diagnostic tool. 

The metamaterial used for this study has been derived 
from the structure previously shown to support a strong 
toroidal dipolar response under linearly polarized excita¬ 
tion HU. Its metamolecule consists of a cluster of four 
rectangular metallic wire loops embedded in a dielectric 
slab [see Fig. [^a)]. The loops are located in two mutu¬ 
ally orthogonal planes, whose intersection gives the axis 
of the metamolecule parallel to ^-axis. Each loop is split 
in the centre of either its top or bottom side depending on 
which of the two orthogonal planes it is located in. The 
metamaterial is formed by translating the metamolecule 
along the x and y axes, resulting in a one-metamolecule 
thick crystal slab with a rectangular unit cell, as shown 
in Fig. [^a). Although the individual metamolecules 
are intrinsically achiral, the resulting metamaterial crys¬ 
tal is chiral due to the regular arrangement of its meta¬ 
molecules (’structural’ chirality m) and is available in 
two enantiomeric forms interconnected by a mirror re- 
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flection [see Fig. 1(a)]. 

The metamaterial sample was fabricated by etching a 
35/im thick copper foil on both sides of a low-loss di¬ 
electric PCB laminate Rogers TMM(R) 3 using high- 
resolution photolithography [see Fig. 1(b)]. The top and 
bottom patterns were electrically connected through nar¬ 
row electroplated holes. All copper tracks were coated 
with a 2/im thick layer of gold to prevent oxidization and 
staining of copper, as well as to reduce Ohmic losses. The 
assembled metamaterial sample consisted of 24 x 24 meta¬ 
molecules. Its transmission response was experimen¬ 
tally characterized at normal incidence in a microwave 
anechoic chamber at frequencies in the range 12.5-16.5 
GHz using broadband linearly polarized horn antennas 
equiped with collimating lenses and a vector network an¬ 
alyzer. Our measurements yielded a complex transmis¬ 
sion matrix tij {i^j = x,^), which fully describes the 
metamaterial transmission in the linear polarization ba¬ 
sis. To obtain the transmission matrix in the circular 
basis we used the following well-known relations: 
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where the subscripts + and — denote right (RCP) and 
left (LCP) circular polarizations. 

As evident from Fig. [^a), the diagonal elements of the 

resulting matrix, and t _, do not match indicating 

the presence of optical activity in the response of the 
metamaterial. The difference between the transmission 

intensities A = — \t _is a measure of structure’s 

circular dichroism and it is particularly strong within two 
resonant bands located at vi = 14.6 and z /2 = 13.6 GHz, 
where the dichroism reaches 60% and 80% respectively 
[see Fig. [2];b)]. 

In order to determine the cause of optical activity at 
these two resonances we conducted a full-wave simula¬ 
tion of the metamaterial response using a 3D finite el¬ 
ement solver (GOMSOL 3.5a). The metamaterial array 
was modelled by placing the toroidal metamolecule in a 
rectangular unit cell with periodic boundary conditions 
imposed along x and y axes [see Fig. 1(a)]. The wires 
of the metallic loops were assumed to be infinitely thin 
strips of perfect electric conductors. Permittivity of the 
dielectric slab was set to e = 3.45 — iO.007, which corre¬ 
sponds to Rogers TMM(R) 3 dielectric laminate used in 
the fabrication of the metamaterial sample. This model 
allowed us to accurately reproduce the main features of 
the metamaterial response, yielding a good agreement be¬ 
tween the calculated and measured far-held transmission 
spectra (solid curves in Fig. [^. The slight discrepancy 
is attributed to fabrication tolerances and uncertainty in 





FIG. 1: (a) A schematic of the toroidal metamolecule and its 
orientatation with respect to the metamaterial crystal lattice 
(indicated by dashed lines). The dimensions are d = 8mm, 
s = 7.5mm, h = 1.5mm, r = 2.44mm, g — w — 0.15mm and 
a = 1.8mm. The metamaterial slab is formed by translating 
the unit cell along the x and y axes, which imposes struc¬ 
tural chirality on the sample. RCP (red) and LCP (blue) 
light propagates along the z-axis of the metamaterial. Mir¬ 
ror image shows the enantiomeric form of the metamaterial 
crystal, (b) Close-up photograph of the fabricated toroidal 
metamolecule, (d) Metamaterial sample after assembly. 


the dielectric constant of the PCB laminate. 

The microscopic origin of the effect was analysed in 
terms of multipolar charge-current excitations supported 
by the metamolecules (analogous to multipolar transi¬ 
tions in conventional molecules). The multipole mo¬ 
ments were computed based on spatial distributions of 
the induced conduction and displacement current densi¬ 
ties extracted from our model, as it was previously de¬ 
scribed in m- To simplify the analysis the expansion 
was truncated at the octopole order for electric and mag¬ 
netic multipoles, and at the quadrupole order for toroidal 
multipoles. Following the recently proposed methodol- 
ogy dS] we can directly link the multipolar excitations 
of individual metamolecules to the far-held response of 
the entire metamaterial array. An adequate accuracy 
of this approach is demonstrated by Figs. [^c) and 
l^d) where we compare the metamaterial transmission 
calculated through the multipole decomposition (dashed 
curves) with that obtained directly from the simulations 
(solid curves). The small quantitative mismatch here is 
the result of limited spatial resolution of extracting the 
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FIG. 2: Measured (left) and calculated (right) circular trans¬ 
mission response of the toroidal metamaterial. Panels (a) and 
(b) show the amplitudes of the direct circular transmission 
and while panels (c) and (d) - circular dichroism 

A. Results calculated using two different methods - finite- 
element simulation and multipole decomposition of currents 
are shown by solid and dashed curves, respectively. Purple 
and grey bands mark the resonances of corresponding con¬ 
ventional and toroidal circular dichroism. 


current data rather than truncation of the multipole se¬ 
ries. 

For the purpose of the analysis it is sufficient to per¬ 
form the multipole decomposition of the metamaterial’s 
polarization conversion response, t^y and tyx^ since these 
can be regarded as a manifestation of the observed cir¬ 
cular dichroism in the linear basis. Indeed, according to 
Eq. (1) the difference between the diagonal elements of 
the circular transmission matrix exists only for non-zero 
linear polarization conversion terms, where \txy\ = \tyx\ 
as dictated by the reciprocity and structural symmetry of 
the metamaterial with respect to the propagation direc¬ 
tion (z-axis). Figures [^a) andj^b) show the results of 
the multipole decomposition where we retained only six 
multipole components for clarity, namely x-components 
of electric p and magnetic m dipoles at z/i, and electric 
quadrupole Qe and ^-components of electric, magnetic 
and toroidal T dipoles at 1^2 (other multipole compo¬ 
nents are at least one order of magnitude smaller than 
the dominant ones): 

p = it 

T=Tfe/[(r-j)r-2r2j]d3r 

ea/3 = 2 b / Ki/3 + (r ' j)] d?r 

where j is the induced current density. 

In agreement with the modelled circular dichroism 
data [Fig. 2(d)], the resonances in the multipole scatter- 
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FIG. 3: Amplitudes of the helds radiated by the four most 
dominant multipole components (electric dipole p, magnetic 
dipole m, toroidal dipole T and electric quadrupole Qe) that 
contribute to polarization conversion response of the meta¬ 
material txy (a) and tyx (b). (c)-(f) give the held distribu¬ 
tions {log\E\ and log\H\) around a metamolecule when ex¬ 
cited by LGP light at ui (c)-(d) and U 2 (e)-(f). The black 
arrows give the corresponding E and H held direction. The 
green (blue) arrows show orientation of the excited electric 
(magnetic) dipoles. 


ing occur at frequencies ui and U 2 . Figs. 3(a) and 3(b) 
indicate a symmetry to the linear conversion response 
though the multipoles contributing to the same resonant 
features in txy and tyx are different. This means that the 
effect of circular dichroism is controlled by essentially 
pairs of different multipolar excitations. In particular, 
at ui the dominant pair comprises electric and magnetic 
dipoles with collinear moments - Px and Such a 

combination corresponds to the textbook optical activity 
and hence is considered as a reference case in our further 
analysis below. At z^ 2 , however, the contributions of Px 
and rUx become only secondary, while the circular dichro¬ 
ism appears to be underpinned by previously unexplored 
and rather exotic combination of the toroidal dipole and 
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electric quadrupole. 

The presence of these multipole pairs can be visually 
detected in the distributions of electric and magnetic 
fields inside the unit cell arising under circularly polar¬ 
ized excitation [see Figs. 3(c) - (f)]. At vi the field lines 
of both electric and magnetic fields are seen to stream 
parallel to x-axis, which clearly correspond to electric 
and magnetic dipolar excitations of the metamolecule 
with the net dipole moments being collinear and oriented 
along x-axis [Figs. 3(c) and (d)]. At z /2 magnetic field is 
seen to confine within a well-defined ring-like area where 
the field lines form a vortex threading all loops of the 
metamolecule [Fig. 3(e)]. Such a magnetic field config¬ 
uration is unique to the toroidal excitation with the net 
dipole moment aligned parallel to the axis of the meta¬ 
molecule (i.e. ^-axis). The distribution of electric field 
features a pattern similar to that in Fig. 3(c) but the field 
lines on the opposite sides of the metamolecule stream 
anti-parallel to each other [Fig. 3(f)], which is typical for 
an electric quadrupole excitation. 

Our conclusions regarding the nature of the observed 
circular dichroism can be further reinforced by examining 
the role of the identified multipole pairs in shaping the 
metamaterial’s polarization eigenstates. The eigenstates 
were computed by solving the characteristic equation of 
the transmission matrix t, which was reconstructed from 
the multipole data. The polarization eigenstates were 
characterised in terms of their ellipticity angle e using 
the following relation: 
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1 

-arcsin 
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where e+ and e_ are eigenvectors of t in circular ba¬ 
sis, and |e| takes values between 7r/4 (circular polariza¬ 
tion) and 0 (linear polarization). A positive ellipticity 
angle defines right handed rotation, while negative - left 
handed rotation. The two limiting cases, e = ±7r/4 
and e = 0, represent polarization eigenstates of optically 
active (chiral) isotropic media and linearly birefringent 
(achiral) anisotropic media, respectively. The calculated 
eigenstates of the metamaterial correspond to the inter¬ 
mediate case of counter rotating elliptical polarizations 
with e ±7r/8 near the resonances, where the ellipticity 
arises due to anisotropy in the metamaterial structure 
[see Fig. |4];a)]. 

The roles of the multipoles are assessed by removing 
their contributions in pairs from the linear conversion re¬ 
sponse, preserving the symmetry of the conversion and 
then recalculating the ellipticity angles of the resulting 
eigenstates. The application of this approach is illus¬ 
trated in Fig. Sb) for the reference case. One can see 
that at vi the removal of just Px and rux leads to a de¬ 
crease of |e| to approximately 7r/60 for both eigenstates 
with their polarizations collapsing to almost linear states. 



Frequency (GHz) 


FIG. 4: Ellipticity angles of the metamaterial eigenstates 
(left) and associated polarization ellipses at vi and 1^2 
(right), calculated based on the multipole expansion of mi¬ 
croscopic charge-current excitations supported by toroidal 
metamolecules. Panel (a) shows the results obtained using 
the complete multipole set. Panels (b) and (c) show the 
results corresponding to pairs px — rux and Ty — Qe being 
removed from the microscopic response of the metamaterial 
(solid curves) and compared to the complete multipole re¬ 
sponse from the upper panel (dashed curves). 


Such a dramatic transition towards purely anisotropic re¬ 
sponse confirms that on the microscopic level the effect 
of circular dichroism at ui results from electric and mag¬ 
netic dipolar charge-current excitations induced in the 
metamolecules. Using the same approach for the combi¬ 
nation of Ty and Qe we also observed a clear change in the 
polarization eigenstates but at z/ 2 , where |e| decreased to 
approximately 7r/30 [see Fig. [^c)]. The somewhat larger 
new values of |e| here compared to the reference case at 
vi are attributed to the non-negligible contribution of Px 
and rrix that also impacts the optical activity at z/ 2 , as 
evident from Fig. Sb). Nevertheless, the resulting pro¬ 
found modification of the eigenstates at z /2 supports our 
initial conclusion that the resonant excitation of toroidal 
dipole and electric quadrupole is the dominant mecha¬ 
nism of the effect of circular dichroism at this frequency. 

The presence of strong toroidal and quadrupole contri¬ 
butions to the optical activity of the metamaterial may 
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seem surprising. Nothing of the kind has been observed 
in chemistry, where the circular dichroism is a standard 
analytical technique that reports, for example, on sec¬ 
ondary and tertiary structures of proteins |20]. Indeed, 
the conventional expression for the rotatory power does 
not feature toroidal or quadrupolar moments m- 

Rnk = Im[{'tpk\P\'>Pn) {'>Pk\m\'tpn)\ , (3) 


reveals that the metamaterial exhibits both conventional 
optical activity, by way of electric and magnetic dipole 
coupling, and, previously unobserved, strong optical ac¬ 
tivity underpinned by the resonant excitation of toroidal 
dipole and electric quadrupole. Such an exotic multi¬ 
pole combination is anticipated to be the cause of (or 
at least strongly contribute to) optical activity in many 
other material systems with toroidal topology, common 
in chemistry and biology. 


and contains only the electric and magnetic dipole 
operators p and ih. This expression is obtained using 
perturbation theory under the common (in chemical sci¬ 
ences) assumption that the size of a molecule is much 
smaller that the wavelength of light (i.e. k • r ^ 1, 
where k is the wavevector). Such assumption, however, 
is no longer valid for the metamolecule, and therefore the 
power series expansion of the complex exponential form 
of the plane wave must be carried out to at least the sec¬ 
ond order of k • r. In particular, for the vector potential 
we have: 


A = AoiJe (4) 

« AoiJe I (^1 + ik-r- | (5) 

Because of the second order term (k • r)^ the magnetic 
perturbation operator A • p [22] acquires the following 
components: 


(k • r)^ V = (k • k)(r • r)V - ([k x r] • [r x k]) V, (6) 
where. 


(r • r)V = T (r 


r X 


r X V 


(7) 


The last term in Eq. [^corresponds to the toroidal mo¬ 
ment. The presence of toroidal contributions in the cir¬ 
cular dichroism is therefore the consequence of non-local 
response of the metamolecules. The complete derivation 
of the expression for the rotatory power that explicitly 
features the toroidal dipole operator T will be published 
elsewhere. 

In summary, we investigated experimentally and theo¬ 
retically a structurally chiral metamaterial that exhibits 
strong circular dichroism and supports non-negligible 
toroidal dipolar response. Analysis of the fields scattered 
by the charge-current multipoles induced in the meta¬ 
molecules and the polarization eigenstates of the system 
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